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Abstract. The emergence of a large scale magnetic field from randomly 
forced isotropic strongly helical flows is discussed in terms of the inverse 
cascade of magnetic helicity and the a-effect. In simulations of such flows 
the maximum field strength exceeds the equipartition field strength for 
large scale separation. However, helicity conservation controls the speed 
at which this final state is reached. In the presence of open boundaries 
magnetic helicity fluxes out of the domain are possible. This reduces the 
timescales of the field growth, but it also tends to reduce the maximum 
attainable field strength. 



1. Introduction 

It was since the mid-seventies when Frisch et al. (1975) and Pouquet et al. 
(1976) came up with the idea that the large scale magnetic fields seen in 
many astrophysical bodies could be caused by an inverse cascade-type phe- 
nomenon. Although there were close links with earlier results that helicity 
and lack of mirror symmetry are important (Steenbeck et al. 1966; see also 
Krause & Radler 1980), the notion of an inverse cascade has put dynamo 
theory into a self-consistent framework of nonlinear turbulence theory. 

Unfortunately the inverse cascade concept was not easily assimilated by 
the astrophysical community. The reasons are simple: the inverse cascade 
approach was developed in the framework of isotropic homogeneous tur- 
bulence, and was not readily applicable to astrophysical bodies that were 
stratified and enclosed in boundaries. Thus, people continued to use a 2 and 
aO-dynamos (Moffatt 1978), which enabled modelling of a large variety of 
astrophysical bodies. 

Here we want is to look more closely at inverse cascade dynamos. In 
particular, we want to know what kind of field they produce and how this 
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relates to the fields generated by an a 2 -dynamo. The full results of this 
work are presented in a separate paper (Brandenburg 2000, hereafter re- 
ferred to as B2000). In the present paper we also discuss the effects of open 
boundaries allowing magnetic helicity fluxes out of the domain into the 
exterior or across the equator. This is an important issue that has been 
raised recently in the context of dynamo theory (Blackman <fe Field 2000, 
Kleeorin et al. 2000). 

In order to model isotropic random flows we adopt a forcing function 
that consists of randomly oriented Beltrami fields with a wavenumber, kf, 
that is larger than the smallest wavenumber, k\, that fits into the box. For 
most of the calculations we use kf = 5 and k\ = 1, leaving some margin for 
scale separation. In order to see more clearly the effects of scale separation 
we also have one run where kf = 30. In the wavenumber band 4.5 < \k\ < 
5.5 (for kf = 5) there are 350 wavevectors which are chosen randomly 
at each timestep, so the forcing is (5-correlated in time, but the resulting 
velocity field is not. In fact, the velocity has a well defined correlation time 
that agrees well with the turnover time r = £f/u rms , where if = 2ir/kf is 
the forcing scale and u rms is the rms velocity. 

The degree of turbulence that develops depends on the range of length 
scales left between the forcing scale and the dissipative cutoff scale. A rea- 
sonable range can only be obtained if the forcing wavenumber is not too 
high, so the run with kf = 30 (Run 6 in B2000) is an example where the 
flow is laminar on the forcing scale. The degree of mixing, as measured 
by the ratio of the turbulent to the microscopic diffusion coefficients for a 
passive scalar, D t /D, is here of order one. In a more turbulent run, Run 
3 of B2000, this ratio is around 40. However, quite independently of how 
turbulent a run is, we find the emergence of a large scale magnetic field 
after some time. The resulting field resembles closely that obtained from 
an a 2 -dynamo with the same (periodic) boundary conditions. This anal- 
ogy enables us to make contact with mean-field theory and to explain the 
resulting turbulent transport coefficients. 

2. Emergence of a large scale field 

The inverse cascade is traditionally described in terms of energy spectra. 
In Figure 1 we compare the spectral field evolution for two different forcing 
wavenumbers. Note that the envelope of the magnetic energy fits under- 
neath a k^ 1 slope. The peaks at k = k\ and kf also fit underneath the same 
slope. There are several features of the spectrum that are characteristic also 
of other cases investigated. For large enough scale separation one sees that 
the magnetic energy grows fast at two distinct wavenumbers, k ss 30 and 
k £3 7. However, when the energy at k ~ 7 reaches saturation the energy 
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begins to be transferred to larger scales until much of the magnetic energy 
is at the largest scale possible. During this phase the magnetic energy at 
intermediate scales decreases to some minimum value which follows roughly 
a /c +3 / 2 spectrum. This effect may be referred to as 'self-cleaning', because 
by removing energy at intermediate scales the field at the largest scales ap- 
pears less perturbed and hence cleaner. This self-cleaning effect is the result 
of nonlinearity, which suppresses the growth at intermediate scales. How- 
ever, the type of nonlinearity does not seem to matter: even with ambipolar 
diffusion the same behaviour is found (Brandenburg & Subramanian 2000). 




Figure 1. Magnetic energy spectra for runs with small and large scale separation, kf — 5 
and 30, respectively. 

In Figure 2 we show cross-sections of B x at different times. Towards the 
end of the evolution the large scale magnetic field is essentially a Beltrami 
field which is here of the form B ~ (siny, 0, cos y), apart from some phase 
shift in y. 

3. The final equilibrium field strength 

Characteristic to all the runs reported in B2000 is the fact that super- 
equipartition field strengths are reached. In Figure 3 the evolution of mag- 
netic and kinetic energies are shown for two runs with different forcing 
wavenumbers. Note also the extremely slow evolution past the moment 



4 



AXEL BRANDENBURG 




-a -2-1 1 3 3 -3-2-1 1 2 3 -3-2-10123 



Figure 2. Gray-scale images of cross-sections of B x (x, y, 0) for Run 3 of B2000 at different 
times showing the gradual build-up of the large scale magnetic field after t = 300. Dark 
(light) corresponds to negative (positive) values. Each image is scaled with respect to its 
min and max values. 



where the kinetic energy drops suddenly to a smaller value. This is when 
saturation of small scale magnetic energy is reached. However, after that 
moment the large scale magnetic energy continues to grow for some time, 
because the resulting large scale field is force-free and does hence not affect 
the velocity field directly. 

The prolonged saturation behaviour found in the present simulations is 
at first glance unusual. Since at late times most of the magnetic energy is 
in the large scales, this slow evolution must have to do with the properties 
of the large scale field. An important property of this large scale field is 
that it possesses magnetic helicity. At the same time magnetic helicity is 
conserved by the nonlinear terms and can hence only change resistively 
on an ohmic timescale. (The case with open boundaries is different and 
will be discussed separately.) In order to demonstrate magnetic helicity 
conservation we consider a periodic box and write the magnetic field, B, 
as the curl of a vector potential, A, so B = V • A. We use the uncurled 
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Figure 3. Evolution of kinetic (dotted) and magnetic (solid) energies. The slow evolution 
of the magnetic energy follows approximately the 1 — e~ 2r,At behaviour (dashed line) of 
Eq. (6) that results from helicity conservation. 



induction equation, 

dA/dt = u x B — rjJ - V</>, (1) 

where J = V x B is the current density, rj is the magnetic diffusivity, and 
4> is the electrostatic potential. The magnetic helicity, which can be defined 
as (A ■ B), satisfies 

d(A- B)/dt = -2r](J ■ B), (2) 

where angular brackets denote volume averages over the periodic domain. 
The steady state solution must satisfy (J ■ B) = 0, so small and large scale 
current helicities must be comparable in magnitude, \(j ■ b)\ f» \(J ■ B}\, 
but of opposite sign. The magnetic helicity is however concentrated on large 
scales (its spectrum is k 2 times that of the current helicity), so its small 
scale contribution is negligible. 

We measure the degree of magnetic helicity of the large scale field by 
the quantity k~^ B = \(A ■ B)\/(B ), which is a length scale characterizing 
the scale of the helical contribution. In a periodic domain of size L the 
smallest wavenumber is ki = 2tt/L, and so k\B is bounded from above by 
kAB < k\ . The large scale current helicity is k\ times the magnetic helicity, 
so we have 

(A ■ B) = T(B 2 )/k AB = (J ■ B)/kl (3) 

where the upper sign applies to the case of positive kinetic helicity in the 
turbulence. Using this in Eq. (2) together with (A ■ B) ss (A ■ B) and 
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(J-B) = (J-B) + (j • b) yields 



Td(B 2 )/dt = ±2r ) k 2 1 (B 2 ) - 2r,k A B{j ■ b). 



(4) 



Prior to saturation (j ■ b) is small, but during saturation its value is limited 
by the kinetic helicity, so 



which is only valid at late times when At = t — i sat > 0. (Here, i sat is the 
time when the small scale field saturates.) We emphasize that this relation 
is rather general and independent of the actual model of field amplification, 
because we used only the concept of magnetic helicity conservation. The 
important point here is that full saturation is only obtained after an ohmic 
diffusion time. In that sense Eq. (6) poses a constraint on the mean magnetic 
field at late times. It applies only as long as the magnetic field is helical. 
Indeed, lower degrees of helicity, i.e. smaller values of k^s, allow larger 
values of the final field strength; see Eq. (6). This has been verified in a 
model where differential rotation or shear contributed significantly to the 
field amplification (Brandenburg et al. 2000). However, this only relaxes 
the constraint by a certain factor which depends on the degree of helicity 
of the large scale field which, in turn, depends on the degree of linkage of 
poloidal and toroidal field. We found that to a good approximation this 
factor is given by the ratio Q of toroidal to poloidal field strengths. In the 
sun this factor is less than a hundred, so this is a relatively minor effect 
compared with the value of the magnetic Reynolds number (10 8 — 10 10 ). 

4. Helicity exchange across the equator 

A different way of relaxing the slow growth problem is to allow for fluxes 
out of the dynamo volume either into the exterior of the dynamo (the 
corona or halo) or from one hemisphere to the other. In any case, there 
would be an extra surface term in Eq. (2). Here we want to discuss the 
latter alternative of a helicity flux between the two hemispheres. Such a 
flux could result from a turbulent exchange of magnetic helicity between 
the two hemispheres and should therefore be proportional to some turbulent 
magnetic diffusivity 7] t . Based on dimensional arguments, one may expect 
such a term to be of the form rjtk 2 s H, where H is the gauge-invariant 



(j ■ b) « (p)(w • u) « ±k i {pu 2 ) = ±k i B'L/n Q - 



(5) 



This yields an evolution equation for the mean magnetic field, 
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magnetic helicity for open volumes (Berger & Field 1984), which replaces 
(A ■ B). The diffusion of magnetic helicity depends on the length scale 
2n/k e f[ over which the magnetic helicity varies (if evaluated over different 
volumes), so we expect k c s < k\. Equation (4) becomes then 

T d(B 2 )/dt = ±27 1 k 2 1 (B 2 ) - 2 V k^ B {j ■ b) ± 2r ]t k 2 cS (B 2 }, (7) 

so the solution for the mean magnetic field is (assuming again (j ■ b) « 
±k f B 2 Jfj, ) given by 

fj - ^$§T {1 " exp [-2 {*l + ( t - W ] } ■ (8) 

Thus, the time dependence is no longer resistively dominated, because the 
microscopic diffusivity is now supplemented by an additional turbulent dif- 
fusivity. Unfortunately, however, the amplitude of the final field decreases 
in such a way that the initial linear growth in unchanged. This is simply 
because of the fact that the flux term, as modelled here, does not act as an 
effective driver, which is what the (j -6)-term did. This is also clearly seen in 
a simulation where we have included an equator by modulating the forcing 
function such that the kinetic helicity varied sinusoidally n the z-direction 
within the domain; see Figure 4. 




Figure 4- Evolution of the magnetic energy for a run with homogeneous forcing func- 
tion (solid line) and a forcing function whose helicity varies sinusoidally throughout the 
domain (dotted line) simulating the effects of an equators. 

5. Conclusions 

The present work has shown that for helical velocity fields a large scale 
magnetic field is generated. There are strong parallels with the fields re- 
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suiting from mean-field a 2 -dynamos (see B2000 for details). One aspect 
that has now begun to receive major attention is related to magnetic he- 
licity conservation, which prevents rapid growth of magnetic helicity and 
hence helical large scale magnetic fields. Shear, which corresponds to dif- 
ferential rotation in a rotating system, relaxes this constraint only partially 
in that it lowers the fraction of the field that contributes to magnetic he- 
licity. On the other hand, by allowing magnetic flux to escape through the 
boundaries, or allowing for mixing of magnetic helicity of opposite sign at 
the equator, the helicity constraint is modified such that the time scale is no 
longer resistively dominated. The problem however is that various attempts 
to model this effect result in significantly lower equilibrium amplitudes of 
the magnetic field. Part of the problem is that the loss of magnetic helicity 
implies at the same time a loss of magnetic energy. It would therefore be 
advantageous for the dynamo to lose preferentially small scale magnetic he- 
licity and energy. This is something that the dynamo in the computer keep 
refusing to do. It may therefore be important to resort to more realistic 
simulations where the flows are driven naturally and not by some artificial 
stirring in space. 
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